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We present azimuthal angle correlations of intermediate transverse momentum (1 — 4 GeV/c) 
hadrons from dijets in Cu+Cu and Au+Au collisions at ^/sJfN = 62.4 and 200 GeV. The away-side 
dijet induced azimuthal correlation is broadened, non-Gaussian, and peaked away from Acp = tt 
in central and semi-central collisions in all the systems. The broadening and peak location are 
found to depend upon the number of participants in the collision, but not on the collision energy or 
beam nuclei. These results are consistent with sound or shock wave models, but pose challenges to 
Cherenkov gluon radiation models. 

PACS numbers: 25.75.Dw 



Heavy ion collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) produce QCD matter at enormous energy 
density [l[ , exceeding that required for a phase transition 
to partonic, rather than hadronic, matter. The produced 
matter exhibits collective motion 0] and is opaque to 
scattered quarks and gluons. The opacity is observed via 
suppression of high momentum hadrons and intermedi- 
ate energy dijets and provides clear evidence of large 
energy loss by partons (quarks or gluons) traversing the 
medium. A key question is how the hot, dense medium 
transports the deposited energy. 

As partons fragment into back-to-back jets of hadrons, 
angular correlations of the hadrons are used to study 
medium effects upon hard scattered parton pairs. 
Hadron pairs from the same parton appear at Acjj ^ 
(the near-side), while those with one hadron from each 
parton in the hard scattered pair appear at A(j) ~ tt (the 
away-side). For brevity, we will refer to these dijet in- 
duced dihadron azimuthal correlations with the abbrevi- 
ation "dijet correlations". 

Of great interest are intermediate transverse momen- 
tum (pt) hadrons, as they can arise from intermediate 
energy jets or involve partons from the medium 0, 0]. 
Their correlations can provide information about energy 
loss mechanisms, dissipation of the radiated energy in the 
medium, and collective modes induced by the deposited 
energy. Theoretical ideas include: Mach cones from den- 
sity waves induced by supersonic partons co-moving 
radiated gluons producing "wakes" in the medium [^, 
ultrarelativistic partons creating Cherenkov gluon radi- 



ation and medium-induced gluon radiation at large 
emission angles 0, S] ■ They all imply significant modi- 
fications of dijet correlations in the away-side, when the 
parton path through the medium is long. In particular, 
some of these theoretical models 0, 0, imply a transi- 
tion from the peaked distribution at Afjj ~ tt character- 
istic of p-l-p and p-t-A collisions to a distribution with a 
peak away from A(f) '--^ tt in head-on Au-|-Au collisions. 

Low pt (> 0.15 GeV/c) hadrons associated with high 
Pt hadrons (> 4 GeV/c) have modified away-side di- 
jet correlations and softened px distributions relative to 
those in p-l-p collisions, suggesting that at least some of 
the lost energy is thermalized in the medium 9]. At in- 
termediate pt, a strong non-Gaussian shape modification 
of the dijet away-side correlation 10] indicates the possi- 
ble existence of a local minimum at A(j> = tt. This Letter 
introduces new parameters to quantify this shape modi- 
fication and reports their dependence on collision energy, 
system size, transverse momentum and centrality mea- 
sured by the PHENIX experiment at RHIC. 

The minimum bias data were collected in the years 
2005 (Cu-I-Cu at ^/IWW = 200 and 62.4 GeV), 2004 
(Au+Au at y^IWN = 200 and 62.4 GeV), and 2003 
(d-f-Au at ^snn = 200 GeV). Charged hadrons are 
tracked using the Drift Chambers and Pad Chambers 
of the PHENIX central arm spectrometers at midrapid- 
ity (hi < 0.35) in the same way as described in [l0|. 
The number of events in the Au-|-Au data at = 
200 GeV used here is about 30 times higher than that in 
[lo| . Collision centrality and the number of participant 
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nucleons (Npart) are determined using the Beam-Beam 
Counters (BBCs) and Zero Degree Calorimeters 

Relative azimuthal distributions Ysamei^4') between 
"trigger" hadrons with 2.5 < px < 4 GeV/c and "associ- 
ated" hadrons with 1 < px < 2.5 GeV/c are formed. We 
correct their shape for the non-uniform azimuthal accep- 
tance of the PHENIX central arms by using the mixed 
event pairs Ymixedi^'P) from the same data sample: 



C(A</>) 



Y, 



3(A'/') jYrr. 

X 



erf(A0)dA(/. 



(1) 



Extensive Monte-Carlo simulations were performed to 
ensure that the true pair distribution shape is recovered 
through this procedure. 

In Au-f-Au and Cu-|-Cu collisions, hadrons have an az- 
imuthal correlation with the reaction plane orientation 
^RP which is proportional to 1 -I- 2u2cos(2(0 — ^rp)). 
This generates a significant correlated background to our 
dijet source J{A(f>) of azimuthal correlations: 

C{A(t>) = 6o(l + 2K''"°'=)(4"^»)cos(2A?i)) + J(A0) (2) 

The charged hadron (U2), where " ()" signifies an event 
average, was measured through a reaction plane analysis 
using the BBCs (3 < \r]\ < 4) as in [13, IH . 

Hadrons have also a much smaller fourth order az- 
imuthal correlation with the reaction plane orientation. 
Its effect was studied with the Au-f-Au data at 200 GeV 
by including the corresponding 2(w4^*°'') {v*"^^) cos(4A0) 
component in the background term of Eq. [2l where the 
(vi) values have also been measured through the reaction- 
plane analysis No significant V4 systematic effects 
on the shape of the dijet correlations were found. 

The background subtraction generates point-by-point 



(A0 dependent) systematic errors from 



trigg\ 



uncertainty and an overall (A0 independent) systematic 
error from 60 uncertainty. The sources of {v2'^^°'^){v*2^^^) 
uncertainty are the (^2) systematic error domi- 
nated by the reaction plane resolution uncertainty, the 
(^2) statistical error, and the systematic error from the 



f^^assoc^ . I^^trigg-^ factorization approxi- 
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mation made in Eq. [2l The latter is estimated to be at 
most 5% of the (W2) product for the most central events. 

The 60 uncertainty is estimated by using three inde- 
pendent methods to calculate 69. The first, called Zero 
Yield At Minimum (ZYAM), assumes that there is a re- 
gion in A0 where the dijet source of particle pairs is neg- 
ligible. 60 is varied until the background component in 
Eq. [5] matches the measured correlation C(A0) at some 
value of A0. In the second method a functional form for 
J(A0) is added to the background, and the sum fitted 
to the measured correlation with 60 as a free parameter. 
Motivated by the theoretical ideas discussed in the in- 
troduction, we use a function that contains a near-side 
Gaussian, and two symmetric away-side Gaussians: 

J( A0) = G( A0) -f G( A0 -71-0) + G{A(j> -n + D) (3) 



While the choice of this functional form is not unique, 
it does provide a reasonable fit to the measured correla- 
tions, as shown by the dotted line in Fig. [TJ The pa- 
rameter D, or peak angle, is motivated by an attempt 
to describe the away-side dijet correlation in terms of its 
symmetry around Aip ^ tt. We note that it also tends to 
absorb any non-Gaussian character of the dijet correla- 
tion. The third method is independent of the measured 

C{A(I)). We calculate bo = ^n{ntrigg){nassoc) / {nsame) 

with hadron production rates measured from all events 
within each centrality class and scale by the same-event 
pair rate, k is a correction for pair-cut bias and ^ is a 
correction for residual correlations due to averaging pro- 
duction rates from events of different multiplicity within 
the same centrality class [l3l |. 

As shown in Table |T] for the Au-I- Au data at 200 GeV, 
there are slight 69 variations depending on which method 
is used to extract its value. However, the resulting shape 
of the dijet correlations is essentially independent of these 
variations. 



TABLE I: bo values in Au+Au data at ./sWn = 200 GeV: 
ZYAM values (first row); variation of fit values from the 
ZYAM values (second row); variation of combinatorial val- 
ues from the ZYAM values (third row). 

Centrality 60-90% 40-60% 20-40% 10-20% 5-10% 0-5% 

ZYAM 60 0.861 0.942 0.960 0.971 0.982 0.988 

fit 5bo -0.003 -0.003 -0.006 -0.028 -0.035 -0.022 

comb. Sbo -0.086 -0.013 -0.004 +0.002 +0.001 +0.001 



Figure [T] summarizes the extraction with the ZYAM 
method of the dijet correlations using the central (0-5%) 
Au+Au data at ^/snn — 200 GeV: the measured corre- 
lation is shown with squares, the background term with 
a full line, and the background subtracted dijet correla- 
tion with circles for values and boxes for the point-by- 
point systematic errors. The systematic errors are cor- 
related since they depend on the same parameter - the 
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){v*2^^^) uncertainty. For clarity, J{A(f)) is shifted 
up by bo, shown with dashed line, hence its amplitude 
should be read from the right axis. We note that, in this 
case, the measured correlation is flat near A(j) ~ tt, even 
before any background subtraction. Due to the cosine 
modulation of the background, a local minimum should 
develop at A^ ~ tt in the dijet away-side correlation. 

Figured] shows a central and a peripheral dijet correla- 
tion for each colliding system and energy. A remarkable 
away-side feature in central and semi-central collisions 
(< 40%) is the peak location away from Acj) = tt, and 
the appearance of a local minimum at A(j) = tt. To quan- 
tify the significance of this minimum in the Au+Au data 
at 200 GeV, we have studied how much (^2*^°'^) (^r*^) 
would need to change for the away-side to be flat. For 
the four most central bins (0-5%, 5-10%, 10-20%, and 20- 
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FIG. 1: (Color online) The measured correlation C{A(j}) 
(squares) and the dijet correlation J{A(ji) (circles with boxes 
for point-to-point systematic errors) in central Au-|-Au colli- 
sions at ^ysNN = 200 GeV. The full line shows the background 
term and the dotted line shows a C( A0) fit with Eqs. ((2)H-([3ll . 
The left axis shows the measured correlation amplitude and 
the right ajcis shows the dijet correlation amplitude. 



40%) it would have to decrease by 85%(5.1cr), 41%(4.2cr), 
20%(2.3cr), and 23%(2.7f7), respectively, where a is the 
uncertainty. 
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FIG. 2: (Color online) (Di)jet correlations (circles with boxes 
for point-to-point systematic errors) in Au-|-Au and Cu-|-Cu 
collisions at ^/snn ~ 62.4 and 200 GeV. Left panels show 
central collisions, while right panels show peripheral collisions. 



is defined here as all Acj) values above the dijet function 
J(A0) minimum, typically one rad. We extract these 
statistics on only the away-side jet peaks in J(A0); possi- 
ble jet-associated flat underlying distributions, which are 
highly sensitive to the uncertainty in bo and precluded 
by the ZYAM assumption, are not included. 

The rms and kurtosis centrality dependence is shown 
in Fig. [3fa). The rms increases with centrality, indicat- 
ing broadening of the away-side dijet correlation, while 
the kurtosis decreases from the value characteristic of 
a Gaussian shape (three), suggesting a flattening of its 
shape beyond an increase in the Gaussian width. 
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FIG. 3: (Color online) Collision centrality, energy, and sys- 
tem size dependence of shape parameters: (a) kurtosis (filled 
symbols) and rms (open symbols); (b) peak angle D. Bars 
show statistical errors, shaded bands systematic errors. 



TABLE II: Dependence of away-side shape parameters on as- 
sociated hadron pt in central (0-20%) Au-|-Au collisions at 
^/sWn = 200 GeV for 3 < Px""" < 5 GeV/c. First error is 
statistical and second error is systematic. 



assoc 

Pt 


D [rad] 


rms [rad] 


kurtosis 


1-1.5 


1.04±0.03±0.03 


1.02±0.02±0.05 


1.68±0.04±0.10 


1.5-2 


1.07±0.04±0.04 


1.06±0.02±0.05 


1.58±0.05±0.10 


2-2.5 


1.05±0.03±0.06 


1.08±0.04±0.08 


1.38±0.11±0.12 


2.5-3 


1.07±0.06±0.06 


1.09±0.07±0.07 


1.35±0.17±0.12 


3-5 


0.88±0.13±0.16 


1.01±0.11±0.14 


1.31±0.23±0.25 



We parameterize the away-side shape change and de- 
viation from a Gaussian distribution by extracting the 
second and fourth central moments around Ac/) ~ tt 
{fin = ((A0 — 7r)"),7i = 2,4), in the standard form of 
the following statistical quantities: the root mean square 
rms = y/jja and the kurtosis = fi^/fjl^- The away-side 



The peak angle D centrality dependence, extracted by 
fitting dijet correlations with Eq. |3l is shown in Fig. 
|3l^b). It is consistent with zero radians in d-|-Au and pe- 
ripheral nuclear collisions, but rapidly grows to a value 
around one radian in central nuclear collisions. Some 
deviation from zero radians of the peak angle may be 
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related to slight non-Gaussian shapes of the dijet corre- 
lations even without medium modification. This can be 
seen in the kurtosis values for d-|-Au and peripheral nu- 
clear collisions, which have values somewhat lower than 
three. For details on dijet correlations in d-|-Au see 14 1. 
The systematic errors in Fig. [3] come exclusively from V2 
uncertainty. No apparent dependence of rms, kurtosis, or 
peak angle D on collision energy or species is observed. 

Tabic HI] shows the dependence of the away-side shape 
parameters on the associated hadron in the Au-|-Au 
data at 200 GeV for a 0-20% centrahty bin, 3 < p^"^'? < 5 
GeV/c, and the following p?^'^"" bins: 1 - 1.5, 1.5 - 2, 
2 - 2.5, 2.5 - 3, and 3 - 5 GeV/c. The peak angle D 
and the rms have no dependence, while the kurtosis 
is consistent with a slow decrease with px- 

Several phcnomenological models for modification of 
the away-side jet have been proposed; all involve a strong 
response of the medium to the traversing jet. Bow shocks 
propagating as sound, or density, waves in the medium 
produce a peak located away from At/) = tt at approx- 
imately the same angle as seen in the data [1, [3]. If 
the peak indeed arises from a sound wave, its location at 
one radian away from the nominal jet direction implies 
a speed of sound intermediate between that expected in 
a hadron gas and quark-gluon plasma A first or- 
der phase transition would cause a region with speed of 
sound identically zero. This region was postulated [1] to 
reflect sound waves and cause a second away-side peak 
located at about 1.4 radians away from A(/) = tt. No clear 
evidence for a distinct peak is seen in our data. 

If the coupling among partons in the medium is strong, 
then the high momentum parton may induce non-sound 
wave collective plasma excitations Q . In the strong cou- 
pling limit the AdS/CFT correspondence was applied to 
calculate the wake of directional emission from a heavy 
quark traversing the medium, where a peak an^le is found 
at values slightly larger than in these data 16]. 

Thepeak may also arise from Chcrenkov gluon radi- 
ation Such a mechanism should disappear for high 
energy gluons, implying that the peak angle D should 
gradually approach zero with increasing momentum of 
associated hadrons. Table HTl shows that this is not sup- 
ported by the data. The medium may induce gluon radia- 
tion at large angles by mechanisms other than Cherenkov 
radiation |7., 8] . Such models can reproduce the observed 
peak if the density of scattering centers is large and the 
gluon splitting sufficiently asymmetric 7]. However, the 
predicted radiation is very sensitive to the treatment of 
geometry, expansion and radiative energy loss framework 
used. Our detailed measurements constrain the options. 

An important issue is whether the density wave cor- 
relations can survive the underlying medium expansion. 
It was shown that the interplay of the longitudinal ex- 
pansion and limited experimental 77 acceptance preserves, 
and even amplifies, the signal of directed collective exci- 



tations [l5|. However, the creation of a shock wave con- 
sistent with our data requires that 75-90% of the jet's lost 
energy be transferred to the collective mode 0, [ll. 



We have presented azimuthal angle correlations of in- 
termediate transverse momentum hadrons from dijets in 
Cu-|-Cu and Au-|-Au collisions at ^/snn = 62.4 and 200 
GeV. The away-side dijet correlation is seen to be broad- 
ened, non-Gaussian and peaked away from A(f) = tt in 
central and semi-central collisions. The away-side shape 
depends on the number of participants in the collision, 
and not on the beam nuclei or energy. The general 
features of the observed shape can be qualitatively ac- 
counted for by a number of phcnomenological models, 
all having in common a strong medium response to the 
energy deposited by the traversing parton. The system- 
atic data presented here provide quantitative tests that 
could discriminate between these models. 
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